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A derivation of the asymptotic expressions for the threshold intensity of laser 
induced thermal scattering in silica microresonator when illuminated with a 
plane wave is present. The calculation of anti-Stokes thermal combination 
frequencies are made for the spherical high Q-factor microresonators. The 
three modes regime of nonlinear interaction is considered. One pump-driven, 
two signal modes, and one mode of temperature relaxation are taken into 
account, satisfying morphology-dependent input and output resonances. 
There are low power thresholds for laser induced thermal scattering at 
morphology dependent resonances, (c) 2008 Optical Society of America 

OCIS codes: 140.3410, 140.6810 

Stimulated thermal Rayleigh and thermal thermal nonlinear effects in silica microres- 

Mandelshtam-Brillouin scattering from bulk onators due to the extremely high Q-factor of 

liquids have recently been reported [1-3,5]. electromagnetic modes coinciding with " Mor- 

The observations have been treated theo- phology Dependent Resonance" conditions 

retically in terms of entropy and localized (MDR) [6, 16, 17]. At the specific threshold 

thermal fluctuations [4]. Rapid thermaliza- intensity of the laser pump and at the res- 

tion of the molecules which absorb at the onant size parameter of the cavity, the sur- 

laser pump, and the subsequent generation face layer of the cavity significantly enhances 

of large density or temperature fluctuation, the internal field inside the cavity at the inci- 

leads to the enhancement of the thermal den- dent resonant wavelength of laser pump, and 

sity fluctuations and increase the anti-Stokes further more efficiently provides optical ther- 

frequency shift in bulk. There are numerous mal feedback from the internally generated 
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wave through optical bistabihty and instabil- 
ity effects [12-15]. The time dependent shift 
of the eigenfrequencies of the "Whispering 
Gallery Modes" (WGM's) can be caused by 
intensity-dependent mechanisms, such as the 
change in the refractive index which is depen- 
dent on laser light intensity and caused by 
an absorption-induced temperature change. 
The time dependent shift of the eigenmode 
frequencies is due to heating by laser power 
absorption, creating periodic thermal oscilla- 
tions (oscillation regime) in spherical high Q- 
factor silica resonators [12, 13] and aperiod- 
ical oscillation (bistability regime) of scatte- 
ring amplitudes, which have been experimen- 
tally investigated in microtoroid resonators 
[14, 15]. These experimental observations of 
the threshold oscillations and gain enhance- 
ments are associated with thermal nonlinear 
MDR and nonlinear interaction of WGM's 
in microspheres. The purpose of this com- 
munication is to report the calculation of 
thresholds of a new stimulated effect, which 
is believed to be attributable to laser in- 
duced stimulated thermal scattering in mi- 
croresonators. The effects of the laser induced 
thermal scattering (LITS) in microresonators 
at MDR conditions has not been discussed 
in the literature. This article reports new 
knowledge relating to LITS, applicable to mi- 
croresonators. The dielectric permeability of 
the resonator substance depends on the tem- 
perature of the resonator and energy in the 
volume of interacting modes which absorbed 
the laser light power. The dielectric perme- 
ability of the microresonator determines the 
relative shifts of the eigenmodes of the mi- 



croresenator. A change in the dielectric per- 
meability leads to a shift in the eigenmodes 
excited by the laser pump in the microresen- 
ator. Furthermore, independently, due to the 
heating of the resonator by the laser pump 
in input MDR conditions, the eigenfrequen- 
cies of the microresenator also have a pe- 
riodic thermal shift. The maximum absorb- 
tion of the laser power and heating occurs 
at MDR frequencies and a positive shift of 
the eigenfrequencies occurs when heat is re- 
leased from the surface of the resonator. Al- 
ternatively, when heat is generated at the sur- 
face of a resonator, the eigenfrequencies of the 
resonator take on a negative shift. The ther- 
mal coupling of the electromagnetic modes 
in a microresonator provides the appearance 
in the scattering spectra of additional anti- 
Stokes thermal frequency shifts. One would 
expect a lower threshold intensity of LITS 
to be caused by thermal mode overlapping 
and MDR conditions simultaneous with Ra- 
man lasing. The lowering of the threshold of 
the Raman laser emission in silica resonators 
(with radius R = 40 fim) has been reported 
[7]. Early theoretical explanations and the es- 
timations of the low power thresholds for Ra- 
man lasing in microspheres have also been re- 
ported [8,9]. Three modes regime of interac- 
tion were considered. The first mode is the 
thermal mode, the second is the signal and 
the third is the pump mode. Maxwell's elec- 
tromagnetic equations and thermal equations 
for one thermal mode were solved by apply- 
ing the methods of slow varied amplitudes for 
the system of ordinary differential oscillation 
equations [10]. If we take into account the 



2 



partial wave amplitudes of the thermal oscil- 
lations, and compute the scattering power of 
a resonator, (introducing Q-factor as the ratio 
of the field energy inside the mode to incident 
power, and multiplying by the leakage rate for 
power threshold) then we set [11]: 

p^^^^^^^^^A/'/fl)' (1) 
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Here: Kijk and Hijk are the thermal 
and electromagnetic mode overlapping coef- 
ficients estimated below, i corresponds to the 
thermal mode, j and k correspond to the elec- 
tromagnetic eigenmodes; Up is the pump fre- 
quency, Qi and Qj are the Q-factors of i and 
j eigenmodes, /x, is the eigenvalue of thermal 
conductivity equation [11], p is the density 
of the resonator, k is the thermal conductiv- 
ity, Cp is the specific heat capacity, R is the 
radius of resonator. The threshold power is 
obtained at MDR conditions for the optimal 
tuning conditions: 2{Tt + t^) = 1 + u'j/ujp, 
where Tt = ki^pi/ R)"^ / pCpUp and Te = l/IQi 
are the relaxation times of thermal and elec- 
trical modes with uof. Electromagnetic decay 
rate of the effective three-mode interaction 
is Te and the thermal relaxation rate is r^. 
The frequency ojf can be calculated by the 
transcendental eigenvalue equations of Mie's 
theory using the resonant size parameter at 
MDR. For the optimal tuning conditions and 
the lowering of the LITS threshold by MDR 
conditions, the thermal anti-Stokes frequency 
is 111]: 

^h-D{^)' (2) 

Here: D = k/pCp. Although two high Q- 
factor circuits are present in the microres- 



onator, the optimal detuning during three- 
mode LITS can be small if the frequency in- 
terval between the resonator modes is approx- 
imately equal to Qt- WGM's efficiently pro- 
vide the thermal feedback for internally ex- 
cited modes at the specific frequency match- 
ing: Up = Us i: ^T- Equations (1),(2) pro- 
vide the information concerning the absorbed 
substance of a resonator as well as the yield 
of the multi-exponential thermal damping of 
any stimulated processes at times compara- 
ble with a thermal oscillation in Eq.(2). The 
mode overlapping coefficients Kij^ and Hij^ 
have a complicated structure, and consist of 
the high oscillating eigenfunctions [11]. The 
mode overlapping coefficients have the follow- 
ing form [10]: H ^ u]a,T, ~ ipCpV)~^ 
and K ~ aT where: a = euj/AirQi, e is 
the dielectric permeability, is the tem- 
perature coefficient of dielectric permeabil- 
ity, \^ is the effective volume of interaction of 
WGM's [16]. The calculations were made by 
using the known values of material parame- 
ters valid for the experiments: density of fused 
silica p = 2.21 [g/cm^], thermal conductivity 
k = 1.4-10^^ [W/cmK], specific heat capacity 
Cp = 0.67 [Ws/gK], a, = 1.45 ■ 10^^ [A'^^], 
D = 9.5 • 10^^ [cm^/s], refraction index- 
n = 1.46, Qi = 10^ [17]. An asymptotic esti- 
mate of the mode overlapping integral is valid 
for a high index of partial wave amplitudes 
of WGM's n ^ X (here x is the resonant 
size parameter) for n surface modes obtained 
from Cauchy-Bunyakowski-Schwarz inequal- 
ity [10]. The mode overlapping coefficients 
K and H were written under the following 
conditions: (1) small losses of WGM's, which 



3 



are appropriate for the high Q-factor of the 
resonator modes; (2) singhng out two inter- 
acting modes within the homogeneous ther- 
mal mode volume frequency interval. One is 
the pump mode and the other is the reso- 
nant signal (anti-Stokes) mode, namely, both 
input and output resonance conditions [6]. 
The input resonance condition is satisfied for 
a broadband thermal detuning of the input 
and output modes, which spans several high- 
Q MDR's, whereas the output resonance con- 
dition is always satisfied, since the bandwidth 
of LITS spans at least several high-Q MDR's. 
Substituting the asymptotic estimation for 
the mode overlapping coefficients K and H 
in Eq.(l) for the threshold power of LITH 
yields: 
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The threshold of LITS is then determined 
by two circumstances: the effective resonant 
heating absorption and the overlapping mode 
coefficients [11]. In Fig.l the dependence of 
the threshold power on the fused silica spher- 
ical resonator radii is shown. The interact- 
ing nondegenerate WGM's "TE",{TEn,n ~ 
x),{x- resonant size parameter) are coupled 
to thermal modes by thermal nonlinearity 
at wavelength of 840 nm. The computed 
threshold input intensities in a silica nonlin- 
ear sphere vary from 2 to 23 fiW for radii 
of 40 to 110 yum, the Q-factor of a "TE" 
resonant mode is Qi = 10^ [17]. As illus- 
trated in Fig. 2, the threshold incident pump 
power Pth of LITS depends on the thermal 
mode order. More quantitatively the thresh- 
old pump power for LITS Eq.(3) is found to 



Fig. 1. The threshold power of STS for fused 
silica microsphere. Radius and thermal mode 
dependence: 1. TE — Tl modes, 2. TE — 
modes, 3. TE — T^ modes, 4. TE — Tj modes. 
The laser pump is mode-locked Ti : AI2O3 
laser operating at 840 nm with the pulse rep- 
etition frequency 82 MHz and the pulse du- 
ration 1 ps. 

be less 50 ^W, 1 < n < 8, for the radius of 40 
fim, well below the threshold of Raman las- 
ing 86 fiW [7] and stimulated Brillouin scat- 
tering in glass sphere, which has the order 
of 160 W [18]. To be specific, the thresh- 
old power of LITS is less then the thresh- 
old of stimulated Raman lasing and stim- 
ulated Brillouin scattering in silica micro- 
spheres, thus the thermal interaction of the 
eigenmodes can be revealed by the spectra of 
the Raman lasing with thermal instability or 
bistability in microresonators. As illustrated 
in Fig.3, the thermal anti-Stokes combina- 
tion frequencies are dependent on the num- 
ber of the interacting modes. The thermal 
anti-Stokes combinational frequency of LITS 
can be varied over a wide range of frequencies 
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Fig. 2. The threshold power of LITS of T„ 
modes with radius of sphere R = 40 fj,m 1. 
TE - T^, 2. TE - T^, 3. TE - T^, 4. TE - 
T^. The pump is at 1.55 nm of external CW- 
diode laser with 300 kHz line- width [7]. 



Fig. 3. The Anti-Stokes combination fre- 
quency of LITS from fused silica microsphere 
irradiated by CW-diode laser pump at 1.55 
fim [7],R = 40 iim 1. TE - T^, 2. TE - T^ 
3. TE - Tl 4. TE - T^ 



from Rayleigh's frequency of shift via thermal 
Rayleigh scattering [3] to the striction com- 
bination frequency of thermal Mandelshtam- 
Brillouin scattering [2,4]. Such a wide range 
of anti-Stokes frequency can be provided by 
systems such ClS db spherical resonator. The 
anti-Stokes frequencies of LITS can occur 
in spectra of Mandelshtam-Brillouin, Raman 
scattering and linear Mie's scattering. The 
frequency can be varied over the wide range 
and depends on the material parameters of 
the resonator as well as both the effective vol- 
ume of WGM and the surface mode overlap- 
ping. In order to have the description of Ra- 
man scattering and lasing in agreement with 
the experiments it would be necessary to take 
into account the thermal modes in the set 
of electromagnetic equations. The difference 
in the threshold power of interacting thermal 
modes enables the separation of these nonlin- 



ear processes from each other. The thermal 
oscillations and the thermal spatial overlap- 
ping of electromagnetic modes are important 
processes for Raman scattering and lasing in 
microresonators. 
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